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Results ol a high statistics study of chromo field distributions around static sources in pure SU(2) gauge theory 
on lattices of volumes 16*, 32'', and 48^ x 64 at /? = 2.5, 2.635, and 2.74 are presented. We establish string 
formation up to physical distances as large as 2 fm. 
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1. INTRODUCTION 

Various effective models originated from 
the type 11 superconductivity scenario of confine- 
ment 1^. Lattice gauge theory, in principle, of- 
fers the laboratory to test such ideas, as it allows 
for ab initio studies from the QCD Lagrangian. 
Quenched calculations have reached the accuracy 
necessary for quantitative studies of the chromo 
field distributions. These are not easily accessi- 
ble, since (a) the energy density carries dimension 
and therefore imposes a lower limit onto the 
lattice spacing a, due to statistical noise, and (b) 
one is forced to work on rather large lattice vol- 
umes to attain source separations sufficient for 
applicability of the string picture. On top of this 
one is faced with the ubiquitous problem of filter- 
ing ground state signals out of an excited state 
background. 

We have carried out a comprehensive investiga- 
tion of flux tube profiles between static qq pairs. 
Exploitation of state-of-the-art lattice techniques 
for statistical noise reduction and ground state 
enhancement enables us to observe string forma- 
tion up to physical distances of 2 fm, correspond- 
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ing to more than 30 lattice sites at f3 — 2.635. 

2. SIMULATION 

We study lattice volumes of 16^ and 32'' at 
P = 2.5 (a w 0.083 fm), 48^ x 64 at /3 = 2.635 
(a w 0.054 fm) and 32" at /3 = 2.74 (a « 
0.041 fm). The scale has been computed from the 
string tension value ^/Ka = 440 MeV. A hybrid of 
Fabricius-Haan heatbath and overrelaxation has 
been applied for the gauge field update. 

The central observables in our present investi- 
gation, the action and energy densities in pres- 
ence of two static quark sources. 
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are calculated by taking the correlation between 
smeared Wilson loops W and plaquettes □, 
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The flux-tube profiles are probed by varying the 
spatial position of the plaquette which acts like 
a chromo electric (magnetic) Hall detector. To 
minimize contaminations from excited states, the 
temporal position is chosen as close as possible to 
the center of the Wilson loop. By averaging two 
(four) adjacent electric (magnetic) plaquettes, we 
arrive at operator insertions that are symmetric 
in respect to a given (spatial) lattice site n. 
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For measurement of the colour field distribu- 
tions we have restricted ourselves to on-axis sep- 
arations of the two sources. All even distances 
R = 2,4,..., R^,^ with i?„,ax = 8, 24, 36 for 
Ls = 16,32,48, respectively, have been realized. 
In order to identify the asymptotic plateau, T 
was varied up to T = 6. The colour field distri- 
butions have been measured up to a transverse 
distance n± = 2 along the entire QQ axis. In 
between the two sources and up to 2 lattice spac- 
ings outside the sources, the transverse distance 
was increased to ,„ax = 6, 10, 15 for the three 
lattice extents Ls — 16, 32, 48, respectively. In 
addition to "on-axis" positions, n = (ni,n2,0), 
we chose plane-diagonal points n = (ni, 712,^2) 
with 712 < 7i±,max/V2- We averaged over various 
coordinates n, exploiting the cylindrical and re- 
flection symmetry of the problem. 
The field measurements have been taken ev- 
ery 100 sweeps at /3 — 2.5, 2.74 and every 
200 sweeps at /3 = 2.635. At {(3,V) = 
(2.5, 16-*), (2.5, 32^), (2.635, 483 x 64), (2.74, 32^), 
8680, 2046, 248, 1480 measurements have been 
taken respectively. 

The temporal parts of the Wilson loops, ap- 
pearing in the colour field correlator, have been 
link integrated while the spatial parts have been 
smeared. The former procedure is applied to im- 
prove the signal to noise ratio, the latter is essen- 
tial for getting early (in T) ground state domi- 
nance. A more detailed description of the lattice 
methods, analysis procedure and results can be 
found in Rcf. [Q. (A data base of colour images 
of the flux tubes and related quantities can be ac- 
cessed via anonymous ftp from wptsO.physik.uni- 
wuppertal.de. The (compressed) .rgb and .ps files 
can be found in the directory pub/colorfiux.) 

3. RESULTS 

Here, we restrict ourselves to some selected top- 
ics. In Fig. |l| the longitudinal action density pro- 
file for a quark separation r sa 1.6 fm is displayed. 
The string, connecting the sources, is (almost) ho- 
mogenious within a region of extent 1 fm. 

The field distributions within the center plane 
of the flux tube have been studied in detail (see 
Fig. H). For small source separations perturbation 
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Figure 1. Longitudinal action density proflle in 
lattice units at /3 = 2.635, r — 30a sa 1.6 fm. 

theory should apply and one might expect the 
distribution to agree with a dipole shape (up to 
lattice artefacts): 
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with S — R/2. n± denotes the distance from the 
qq axis. Fits to the above parametrization have 
been performed, in which 6 has been treated as a 
free parameter. At large (compared to the width 
of the flux tube) source separation, the string 
picture should apply and we expect a Gaussian 
shape, at least for small n^: 
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Due to cancellations of electric and magnetic 
components, precise data for the energy density 
have only been obtained for r < 0.75 fm. Within 
this region, where the string is not yet developed 
the data is well described by leading order lattice 
perturbation theory with an effective coupling, 
increasing with R in agreement with asymptotic 
freedom. The action density data allow for fits 
up to distances of 2 fm. Ansatze eqn. (^) and (^ 
describe the data equally well as can be seen from 
Fig. I 

The root mean square width, (5, of the flux 
tube is physically more interesting than the half 
width which can be read off from plots like Fig. 



3 



0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 


-0.0005 



p 


R=12 

T (3) — - 
\ (4) 


-fif--- 


:'|if.... 





-10 







10 



Figure 2. Central transverse action density profile 
in lattice units at /3 = 2.5, r = 12a « 1.0 fm, to- 
gether with dipole (3) and Gaussian (4) fit curves. 



In principle, 5 can be extracted from fits but 
strongly depends on the underlying parametriza- 
tion. As we are mainly interested in the R de- 
pendence, especially in probing the predicted log- 
arithmic asymptotic increase of 5{R) [|l|, we use 
a geometric definition of the width [Q that avoids 
large errors on ratios of (5s, obtained at different 
quark separations. It is based on the assumption 
that the transverse shape of the distribution is 
independent of R which holds true to the right 
of the vertical line in Fig. ||. The model depen- 
dence of 5^ has been cast into a (multiplicative) 
geometry factor 7 which is of order one. 

The width increases rapidly until it saturates 
at r « 1 fm at a value 0.5 fm< 5a < 0.75 fm 
where the error is mainly due to the uncertainty 
in 7. Combining the string expectation, 6^ = 
5q ln(i?/i?o)i with our data, yields an inverse cut- 



off wavelength R^' « SOVK (full curve) when 
i5o is constrained to the universal value So = 
l/{TrVK) 0. Without restrictions on 60 we find 
the limits 60 < 2.2 /{nVK) and R^^ > SVK, re- 
spectively (dashed curve) . In case of the simpler 
case of 3d Z2 gauge theory, Gliozzi has observed a 
more pronounced logarithmic increase of S{R) 



4. CONCLUSIONS 

We have demonstrated that Wilson loop- 
plaquette correlations offer a viable access to a 
lattice study of the fiux tube problem on the re- 
quired length scale of one to two fm. The crucial 
ingredient of our method is the smearing of the 
parallel transporter within the bilocal QQ cre- 
ation operator. This secures a controlled ground 
state preparation of long fiux tubes within few 
lattice time slices. 

As a result, we can observe fiux formation of 
the action density over source separations as large 
as 2 fm on lattices with resolution .05 fm and ex- 
tent 2.7 fm. A logarithmic increase of the width 
as suggested by string pictures for the "asymp- 
totic" R region is consistent with our data, sug- 
gesting a surprisingly large ultra violet cut-off 
on inverse wavelengths in effective string models, 
> 1 GeV. 




Figure 3. Width of the action density flux tube, 
S, against R in units of the string tension. 
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